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Abstract — The construction of codes in the projective space for 
error control in random networks has been the focus of recent 
research. The Pliicker coordinate description of subspaces has 
been discussed in the context of constant dimension codes, as 
well as the Schubert cell description of certain code parameters. 
In this paper we use this classical tool to reformulate some 
standard constructions of constant dimension codes and give a 
unified framework. We present a general method of constructing 
non-constant dimension subspace codes with respect to minimum 
subspace distance or minimum injection distance as the union 
of constant dimension subspace codes restricted to selected 
Schubert cells. The selection of these Schubert cells is based 
on the subset distance of tuples corresponding to the Pliicker 
coordinate matrices associated with the subspaces contained in 
the respective Schubert cells. In this context, we show that a 
recent construction of non-constant dimension Ferrers-diagram 
rank-metric subspace codes is subsumed in our framework. 

I. Introduction 

A linear multicast network coding model is defined as 
follows: A source node transmits n packets, each an m-ary 
symbol over F g ; each node in the network transmits F 9 -linear 
combinations of the incoming packets at that node, and at any 
destination node the received packets may be represented as 
the rows of an iV x m matrix Y = AX + DZ, where the 
matrices X, Z are of respective dimensions n x m,t x m over 
¥ q . The matrix X has n source packets as its rows, A is the 
transfer matrix of the network, the rows of Z are the error 
packets, which may be random or introduced by an adversary, 
and D the transfer matrix with respect to these error packets. 
In non-coherent or random network coding, it is generally 
assumed that the source and the destination nodes have no 
knowledge of the network topology - the only assumption 
is that the number of error packets, i.e. the parameter t, is 
bounded. 

The motivation for construction of subspace codes for random 
networks over the projective space came from the seminal 
work of Kotter and Kschischang [1], where the problem of 
error-and-erasure correction in random linear networks was 
translated into that of transmission and recovery of linear 
subspaces of a projective space. Subsequent research (0, 0) 
constructed subspace codes, with subspace distance (0) as 
the metric, based on rank-metric codes. A very significant de- 
velopment in the construction of constant dimension subspace 
codes was the use of Ferrers diagram representation of rank- 
metric codes which fit into a row-reduced echelon form rep- 
resentation of subspaces in a Grassmannian (0), which was 



further employed in conjunction with a lexicographic ordering 
of subspaces ([4]). The above technique and its variants ([6|) 
have resulted in some of the best-known constant dimension 
subspace codes till date. Other important constructions of 
constant dimension subspace codes include those presented 
in 0, 0, @. 

Restriction of codeword subspaces to a particular Grassman- 
nian yields constant dimension subspace codes; a natural 
generalization was to attempt the construction of non-constant 
dimension subspace codes over the entire projective space. In 
this context, a new metric, namely the injection distance, was 
defined in [10], which was shown in some cases (for instance, 
a worst-case adversarial model) to give a more precise measure 
of the performance of a non-constant dimension code than 
the subspace distance metric. In ifTTI . subspace codes (non- 
constant dimension) were constructed for the injection metric 
based on the methods in |3|. Bounds on the size of projective 
space codes, analogous to classical coding theoretic ones, and 
some explicit constructions of non-constant dimension codes 
were presented in lfl2l . 

A new paradigm in random network coding was introduced by 
the so-called orbit codes iTPJl . which generalizes some of the 
results in [7| and describes constant dimension subspace codes 
as orbits of the action of suitable subgroups of the general 
linear group on a Grassmannian. Subsequent work ([14], [15], 
|fl6ll ) characterized orbit codes which result from the action 
of cyclic subgroups. Recently, in [17|, the classical tool of 
Pliicker coordinates was employed to compute the so-called 
'Pliicker orbits' associated with an orbit code and the defining 
equations of Schubert varieties were used in the description 
of constant dimension subspace codes. 
The contributions of this paper are listed below. 

(1) We have used some classical results (|18|) to provide a 
unified framework for some standard constructions of sub- 
space codes (e.g. in 0, 0, 0) in terms of the Pliicker 
coordinates of the constructed subspaces. In particular, we 
show that a strong correspondence exists: the generat- 
ing matrices in both the so-called 'lifting construction' 
and its generalization correspond exactly to the matrix 
of the independent Pliicker coordinates with appropriate 
constraints on the choice of indices. 

(2) We show that the characteristic representation of a Schu- 
bert cell in a projective space in terms of the indices 



corresponding to a non-zero Plucker coordinate can be 
used as an "identifying vector"(|3], j4|, ifTTl ) for sub- 
spaces in that Schubert cell. In particular, we show that 
both the subspace distance and injection distance among 
the subspaces in a projective space can be bounded using 
the subset distance among the defining index sets of 
corresponding Schubert cells. 

(3) We give a construction of non-constant dimension sub- 
space codes which is an extension of a construction 
in |[T3l to the non-constant dimension case - for both 
subspace distance and injection distance. The codes are 
described as unions of Schubert cells selected to guarantee 
a minimum subspace or injection distance. Moreover, we 
show that the Khaleghi-Kschischang construction fTTTl 
for non-constant dimension Ferrers-diagram rank-metric 
subspace codes is subsumed in our proposed construction 
if the constant dimension subcode within each Schubert 
cell is a suitable rank-metric code. 

The organisation of this paper is as follows: Section |TT| is a re- 
vision of pertinent mathematical results regarding the Plucker 
coordinate description of projective subspaces and Schubert 
cells. In Section |IIIJ we reformulate standard subspace code 
constructions in the language of Plucker coordinates and give 
a unified framework. The next section (Section UVb gives 
the general construction of non-constant dimension subspace 
codes. We conclude with a summary of the results and a 
discussion on the directions of ongoing and future research 
in this context. 

II. Preliminary Concepts 

In this section we briefly review and organize some math- 
ematical concepts which will be used in interpreting recent 
results in random network code constructions and lead to the 
formulation of new coding schemes. 

A. Projective Space over a Field 

We begin with the definition of a projective space over an 
arbitrary field. 

Definition 2.1: [19| The projective space associated with a 
vector space V over a field K is the set P(V) of lines in V. 
Moreover, the projective space P(V) is associated with the 
if-vector space V in the following canonical way |fl9l : 
Consider the equivalence relation ~ on the set V* := V^\{0}, 
where, for x, y G V* , x ~ y if and only if y = ax. a G 
K, a ^ 0. Then the canonical map V : V* i — ► P(V) 
associates each vector x with the projective point Kx. 
The dimension of P(V") is defined as Dim(P(F)) := 
Dim(y) — 1; P(0) is empty, and its dimension is — 1. 
Therefore, a d-dimensional subspace or a d-space of the n- 
dimensional projective space P" (K) over a field if is a set of 
points whose representative vectors alongwith the zero vector, 
form a d + 1-dimensional subspace of the n + 1-dimensional 
if -vector space V, where P n (K) = V(V) lEOl . 



B. The Plucker Coordinates of Subspaces 

A d-dimensional linear space Ld in P™ is a set of points 
P = (p(0), • • • ,p(n)) G P" whose coordinates p(j) satisfy a 
system of n — d independent linear equations: 

n 

^2 a ij P(j) = 0; i = 1, • • • , n - d. 

i=o 

There exist d + 1 points Pi — (Pi(0),--- ,Pi(n)) ;i = 

0, ■ • • , d, which span the d-space Ld- Consider the (d + 
1) x (n + 1) matrix P = denote as D(j ,-- - ,j d ) 
the determinant of the (d + 1) x (d + 1) submatrix formed 
by those columns of P indexed by the sequence of integers 
jo < • • • < jd, < jk < 7i. It follows that at least one of 
the N + 1 such determinants, where N := (^i) — 1> must 
be non-zero. Written in a lexicographic order on the column 
indices, (• • • , D(jo, • • • , jd), • • • ) determines a point in P N ; 
these determinants constitute the Plucker coordinates of the 
d-dimensional linear (projective) subspace Ld- 

Theorem 2.2: ([18|, |21|) There exists a natural bijective 
correspondence between the d-spaces in P" and the points of 
P N whose Plucker coordinates satisfy the following relations: 

d+l 

^(-l) Q D(jo,--- ,j d - 1 ,k a )B(k ,--- ,k a ,-- - ,kd+i) = 0. 

Q=0 

where k a denotes omission of k a and < jp, fc 7 < n. 

It is easily verified that the Plucker coordinates of any d-space 

in P™ indeed satisfy the above relations. Proving the converse, 

1. e. any point of P N whose coordinates satisfy the above 
relations corresponds to a unique d-plane in P™, offers some 
interesting insights [ 18"]. For instance, it can be shown that 
any d-subspace of the projective space P™ has an associated 
matrix representation having the (d+l) x (d+l) identity 
matrix as a sub-matrix. The following proposition spells this 
out in detail. 

Proposition 1: [18] There is a natural bijective correspon- 
dence between the set of points of P N whose coordinates 
satisfy the quadratic relations in Theorem 12.21 with the condi- 
tion D(ko, ■ ■ ■ , kd) 7^ 0, and the affine (d+l)(n— d) space of 
(d+l) x (n+1) matrices (pi(j)), i = 0, ••• , d, j = 0, ••• , n, 
such that the (d + 1) x (d + 1) sub-matrix Pi(k a ), i = 
0, ••• ,d, q = 0, ••• ,d is the identity. Moreover, such a 
matrix Pi(j) corresponds to the point of P N with coordinates 
D(ioj • • • ,jd) — Det(pi(jp)); such a point of P N corresponds 
to a (d + 1) x (n + 1) matrix with entries: 

Piti) = D ( fc o>" - ,h-i,j,k i+1 ,--- , fe d )/D(fc , • • ■ ,kd). 

(2) 

It follows from the proposition that the set of points of P N , 
whose co-ordinates satisfy the quadratic relations, is covered 
by N + 1 copies of the affine (d + l)(n — d) space; this 
is the Grassmann manifold or the Grassmannian denoted by 
Q(d + 1, n + 1). In the next section we show that the above 
proposition is the key to the interpretation of some well- 
known constructions for constant dimension codes for random 
networks in terms of the Plucker coordinates of subspaces. 



Henceforth we term the matrix (pi(j)) as the Pliicker coordi- 
nate matrix of the subspace whose representation as a point 
of P N is given by: (■ • • ,Pi(j),- ■ • ) in lexicographic order. 

C. Schubert Cells and their Representation 

Let Q(d,n) be the Grassmannian of d-dimensional sub- 
spaces of an n-dimensional K -vector space. For U £ G(d,n), 
there exists a uniquely determined basis {vi, v%, • • • , Vd}, Vi £ 
K n , of the form: 
Vi = (*,••• ,*,1,0,--' ,0), 
v 2 = (*,■■■ ,*,(),*,••• ,*,1,0,--- ,0), 

Vd = (*,••• ,*,0,V ,*,0,*, ••• ,*,1,0, ••• ,0). 

S v ' 

(d— 1) zeros 

The subspace [/ is uniquely determined by the placement 
of the terminal l's in the Wj's, indexed by a sequence of 
integers 1 < b\ < 62 < • ■ • < bd < n, and the elements 
of K which occupy the * positions. All the subspaces which 
are parametrized by a particular tuple (3 = (&i, 62, • • • ,bd), 
denoting the positions of the trailing l's in a basis of the above 
form, are said to constitute the Schubert Cell Sp of dimension: 

d P = (b 1 -i) + --- + (b d -d) = j2b l -( d+1 ). 

i=l ^ ' 

An alternative description of Schubert cells involves subspaces 
whose basis vectors form a d x n matrix in row reduced 
echelon form, with the defining d-tuple indexing the positions 
of the leading l's instead. The Schubert cell constituted of 
such subspaces, say, S a with a — (ax, 02, • ■ ■ , ad) and the a% 
satisfy 1 < a\ < a 2 < ■ ■ ■ < <Zd < n, has dimension: 

d a = (n — ai — (d — 1)) H h (n — a d ) 




In both the cases, the maximum dimension n(n—d) is achieved 
for the tuple (1, • • • , d) and the minimum dimension is for 
the tuple (n — d + 1, • • • , n). The second characterization of 
Schubert cells is relevant to the construction of non-constant 
dimension subspace codes for random networks. 

III. A Unified Framework for Lifting 
Constructions using Plucker Coordinates 

In this section, we reformulate two well-known methods 
for constructing constant dimensional subspace codes in terms 
of the Plucker coordinates of certain subspaces, providing a 
unified framework for both of them. 

A. Overview of the Lifting Constructions 

The Lifting Construction was first proposed in [ 1 1 to yield 
constant dimensional Reed-Solomon-like codes for error and 
erasure correction in random networks. In a subsequent work 
J21, constant dimension subspace codes were constructed by 
'lifting' rank-metric codes as follows. 



Definition 3.1: [2| Let X be an I x m matrix over ¥ q . We 
define the 'lifting' map as follows: 

I : F^ xm — ► Q(l, l + m); X^ I(X) = ([I X}) 

where / denotes the I x I identity matrix and (M) denotes 
the row-space of the matrix M. The subspace l(X) is called 
a 'lifting' of the matrix X. 

If C is a rank-metric code (e.g. [23 1), the subspace code 1(C) = 
{1(C) I C e C} is called the 'lifting' of C. 
A generalization of the lifting construction was achieved in 
and ifTTI . where matrices in the row reduced echelon form 
(RREF) have been used in conjunction with Ferrers diagram 
representation of rank-metric codes to construct the desired 
subspaces. We now briefly state this construction through some 
relevant definitions fll3l. IfTTI '). 

Definition 3.2: A Ferrers diagram E2ll is a graphical 
representation of integer partitions as an array of dots such 
that the i-th row has exactly the same number of dots as the 
i-ih element in the partition. An mxn Ferrers diagram has m 
dots in the rightmost column and n dots in the topmost row. 

Definition 3.3: Let T be a / x m Ferrers diagram. A rank- 
metric code C is called a Ferrers diagram rank-metric code 
(FDRM code) associated with T if all the codewords are I x m 
matrices in which the positions of the non-zero entries match 
those of the dots in T. In this case, the codeword matrices are 
said to fit the Ferrers diagram T. 

Definition 3.4: Let v £ {0, 1}™ +1 be a vector of weight k. 
The profile matrix of v, denoted P„, is a k x (n + 1) matrix 
in RREF such that: 

(i) The leading coefficients of the rows of P„ appear in the 
columns indexed by supp(w). 

(ii) The entries in P„ other than 0's and l's are dots. 

Extracting the dots of P^ forms a Ferrers diagram which is 
denoted by J-(v). 

Definition 3.5: Let v £ {0, l} n+1 be a vector of weight 
k and let X £ ¥ l q xm , I < k, m < n — k + 1, be a matrix 
which fits J-(v). Then the (generalized) lifting of X (by v) is 
defined as I V (X) = (P V (X)) C G(IJ + m), where P V {X) is 
the profile matrix of v with the dots replaced by the entries 
of X. 

As before, if C is a rank-metric code, then l v (C) = 
{l v (C) I C £ C} is called the generalized 'lifting' of C, which 
is an FDRM subspace code associated with T(v). 

B. Unified Framework Using Plucker Coordinates 

To interpret the above constructions of constant dimensional 
subspace codes in terms of the Plucker coordinates of the 
associated subspaces, we prove the following proposition as a 
consequence of Theorem 12.21 and Proposition [TJ 

Proposition 2: If X £ f( d+1 ) x ( n+1 ' j s j n row-reduced 
echelon form (RREF), then the Plucker coordinate matrix 
associated with the subspace (X) C Q(d+ 1, n+ 1) coincides 
with X, upto multiplication by a unit. 

Proof: Suppose the columns ko, ki, ■ ■ ■ , kd of X contain the 
leading elements of the rows. As D(feoi"" i^d) 7^ 0, by 



Proposition Q] (X) is identified with a point of P N , where 
N = (^i) — 1) corresponding to a (d + 1) x (n + 1) matrix 
with entries: 



is: 



Mi) = D(fco,-- - ,ki-i,j,k 



i+l> 



, fcd)/D(fc , • • • j fed), 



which is the associated Pliicker coordinate matrix. ■ 
Now we give the following interpretations of the lifting 
construction and its generalization in terms of the Pliicker 
coordinates. 

Lifting Construction: In the setting of Proposition [2] set 
d = I — l,n = I + m — 1 and < ko < ■ ■ ■ < kd < n, 
with D(fco,-- - ,k d ) ^ 0. Let X 6 F^ xm , for instance, a 
codeword matrix of an I x m rank-metric code over ¥ q . Select a 
subspace of Q{d+ 1, n + 1) whose Pliicker coordinates satisfy 
the following conditions: fcj = = 0, • • • , d, i.e. we have 
D(0, ■■■ ,1 — 1) = 1 and all other coordinates are obtained 
from those of the form D(j , • • • ,3d)) < j/3 < n, at most 
one of the jp ^ {0, • • • , I — 1}, and that particular column 
is from the matrix X in ascending order. This subspace is 
precisely the lifted subspace l(X) = ([I X]}. 
Generalized Lifting Construction: The generalized lifting con- 
struction may similarly be described as the selection of sub- 
spaces with particular constraints on their Pliicker coordinates. 
Given a vector v as described, set k = d + 1 and let 
hi, i = 0, ■ • • , d, < fcj < n, denote the positions of non-zero 
components of v in ascending order. If X S ¥ l q x m be a matrix 
which fits J-(v), where I < k, m < n — fc + 1, (e.g. X is a 
codeword of a FDRM code associated with the Ferrers diagram 
■F(v)) construct the augmented matrix A(X) G ^kxn-k+i m 
two steps as follows: 

(1) To each column of X append k — I zeros from below to 
form X' G ¥ k q xm . 

(2) Add n — k — m + 1 all-zero columns among the columns 
of X 1 preserving the order of those columns in P v to form 

A(X). 

Select a subspace of Q(d+1, n+1) whose Pliicker coordinates 
satisfy the following constraints: D(fco, • • • , kd) = 1 and the 
remaining coordinates are obtained from those of the form 
D(io7 • • • , jd), < jp < n, where at most one of the 
j/3 ^ {^0) ■ • • > k d }, and that particular column is from the 
augmented matrix A(X) in ascending order. This subspace is 
the lifted subspace l v (X) = (P V (X)). 

We now give a simple example illustrating the Pliicker coor- 
dinate description of the generalized lifting construction. 
Example Let v = 110101 e {0, l} 6 , i.e. n = 5, k = 4. Then 
the profile matrix of v and the associated Ferrers diagram are 
given by: 



1 • • 

1 • • 

1 .0 

1 

In this case, I = 3 < fc = 4, but m = 2 = n — fc + 1. Now if 
X e ¥l x2 is a matrix which 'fits' T(v), then the lifted matrix 



lv(X) = 



10 10 10 

110 10 

1 1 

1 



The positions of the non-zero entries in v are indexed as: fco = 
0, k\ = 1, &2 = 3, &3 = 5. Hence the matrix corresponding to 
the Pliicker coordinate D(0, 1, 3, 5) is the 4x4 identity matrix 
and so, D(0, 1, 3, 5) = 1. As k — I = 1 and m = n — k + 1, 
both the matrix X' and the augmented matrix A(X) are given 
by: 



X' = A{X) = 



Consider the four coordinates: D(2, 1, 3, 5), D(0, 2, 3, 5), 
D(0, 1, 2, 5), D(0, 1, 3, 2), involving the first column of A(X) 
arising out of the column of with index 2. We have: 
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= 1, Det 
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which is in accordance with Proposition |2j and the Pliicker 
coordinate matrix coincides with the lifted matrix I V (X). 
It follows that both the lifting construction and its general- 
ization have precise interpretations in terms of the Pliicker 
coordinate description of subspaces. In fact this interpretation 
provides a unified framework as the encoding procedure is 
shown to be equivalent to a choice of a set of column indices 
from the underlying (codeword) matrices and computing the 
Pliicker coordinates for those particular indices. In the lifting 
construction the 'pivotal' set of indices {ko,-- - ,kd} was 
constrained to be just the set {0, • • • , d}. Hence it actually 
limits the choice of subspaces to the so-called principal 
Schubert cell, with a corresponding matrix representation 
having an identity matrix in the leftmost d + 1 columns. In 
the generalized lifting construction, the (d + 1) x (d + 1) 
identity sub-matrix is interspersed with columns from the 
underlying (codeword) matrix; hence, other Schubert cells 
are also exploited. But, in both cases the lifted matrix still 
coincides with the (d + 1) X (n + 1) matrix formed by the 
'independent' Pliicker coordinates in the setting of Theorem 
12.21 and Proposition Q] as has been outlined earlier. 

IV. A Construction of Non-constant Dimensional 
Subspace Codes 

The construction of non-constant dimensional subspace 
codes has been addressed in and ifTTI : in the former, 
by 'puncturing' a constant dimension code that results from 



a multilevel code construction; in the latter, by partitioning 
the projective space followed by a two-level selection of 
subspaces. In this section we present a construction of non- 
constant dimension subspace codes which is an extension of 
the construction of orbit codes for random networks [fl"3l for 
constant dimension codes. The Khaleghi-Kschischang con- 
struction 1111 also employs the so-called injection distance 
metric, rather than the subspace distance metric, as the former 
sometimes outperforms the latter in a worst-case adversarial 
model. So we first outline our construction for subspace 
distance and then discuss the modifications required with 
injection distance as the metric. 

V. Codes with Subspace Distance as Metric 

First we describe the selection of subspaces in terms of 
their Pliicker coordinates. Recall from Proposition Q~]that every 
d-space in the projective n-space has a non-zero Pliicker 
coordinate D(ko, ■ ■ ■ , kd) such that the (d + 1) x (d + 1) 
submatrix of the coordinate matrix formed by the columns 
indexed by hi, i = 0, • • • , d, is the identity. Let U, V be two 
different subspaces, of possibly different dimensions d u , d v , 
with the tuples {k u } and {l v } as the column indices of the 
identity submatrix in each case. Then the symmetric distance 
between the tuples {k u } and {l v } gives a lower bound on the 
subspace distance d s (U, V) as follows. 

Proposition 3: Let U, V be subspaces of P n of projective 
dimensions s and t respectively; and let {fco, fci, • - • , k s }, < 
k a < n, and {lo,h,--- Jt}, < lp < n, be the tuples 
specifying the identity submatrices of the respective Pliicker 
coordinate matrices. Then d s (U, V) > A({k a },{lp}), where 
A({u},{v}) is the symmetric distance between the sets 

{«}, {«}• 

Proof: Without loss of generality assume that the Pliicker 
coordinate matrices M(U), M(V), of respective dimensions 
(s + 1) x (n + 1) and (t+ 1) X (n+ 1), are in RREF, keeping 
the columns of the identity submatrices intact. Now 

A({fc Q }, {lp}) = \{k a } \ {l p }\ + \{l p } \ {k a }\. 

Let a := \{k a } \ {lp}\, b :— \ {k a }\. Then there are 

s + 1 — a = t + 1 — b matching positions in the matrices 
M(U), M(V) where there are leading l's; hence we have 
Dim(C/ HV) < (l/2)(s +t-(a + b)). But then 

d s (U,V) := Dim U + Dim V - 2Dim([7 n V) 
> s + t-2((l/2)(s + t-{a + b))) 
= a + b. 

■ 

The above result is a general reformulation of Lemma 2 in (3), 
which relates the subspace distance between the subspaces of 
a projective space with the Hamming distance between their 
binary profile vectors. 

Following ifTTl . a two-step strategy for constructing non- 
constant dimension subspace codes over a projective space 
is outlined below: 



(1) Selection of Schubert cells whose constituent subspaces 
across the cells are at a prescribed minimum subspace 
distance from each other; 

(2) Selection of subsets of subspaces in each Schubert cell 
which have the prescribed minimum subspace distance 
among them. 

In IfTTl . the selection of subsets within a Schubert cell is 
achieved via lifting of Ferrers diagram rank-metric (FDRM) 
codes, while a profile vector selection algorithm is imple- 
mented to choose appropriate Schubert cells which also sup- 
port the largest lifted FDRM codes. As discussed in preceding 
sections, Schubert cells can be characterized by the tuples 
which index the positions of leading 1 's of a matrix in RREF. 
Following Proposition [3] the selection of a Schubert cell is 
translated into the selection of a tuple of column indices in 
the Pliicker coordinate matrix of a representative subspace. 
For the selection of subspaces within a Schubert cell we 
propose a method exploiting the characteristic representation 
of subspaces in the Schubert cell. First we state some relevant 
definitions. 

Definition 5.1: Let the set of basis vectors of any subspace 
of a Schubert cell S a have a d x n matrix representation in 
RREF with the positions of the leading l's indexed by the 
tuple {a} = {<2i,a2, • • • ,dd}- The Schubert cell matrix A a 
corresponding to the Schubert cell S a is the d x n matrix in 
RREF in which the columns indexed by the tuple {a} form a 
d x d identity matrix and all other columns are zero. 

Definition 5.2: Given a Schubert cell S a characterized by 
the tuple {a} — {ai,a2, ■•• ,ad}, 1 < a% < ■ ■ ■ < 
o-d < n, the complementary Schubert cell S^, is characterized 
by the tuple {(3} = {b u ■ ■ ■ , b (n _ d) } = {l,2,---,n}\ 
{ai,a 2 , • • • ,a d }. 

Evidently the Schubert cell matrix for is given by the 
(n — d) x n matrix A c a in RREF with the columns indexed 
by {(3} forming an (n — d) x (n — d) identity matrix and the 
remaining columns all zero. 

Let S a be a Schubert cell with the associated d x n Schubert 
cell matrix A a . Also denote the d x n RREF 'matrix' rep- 
resentation of the S a with asterisks indicating arbitrary field 
elements (vide Subsection IH-O as M(S a ). 
Consider the set of d x (n — d) matrices such that the 
columns of these matrices correspond to those columns in- 
dexed by the set: (3 = {bx, b 2 , ■•■ , b n - d } := {1, 2, • • • ,n}\ 
{cti, 02, ■ • • , a d } in ascending order, having non-zero elements 
only in the positions of the asterisks in M(S a ). Select an 
additive group Q of these matrices such that each element of 
the group has rank at least r/2. 

For each G G Q construct the d x n matrix C(G) as follows: 

C(G) — A a + G ■ A c a . 

The constant dimension subspace code constructed within the 
Schubert cell S a is then given by: 

C a = {(C(G))\Geg}. 

The minimum subspace distance of the code satisfies d m i n > 
r, as borne out in the following theorem. 



Theorem 5.3: Let G\, G2 belong to an additive group of 
d x (n — d) matrices with minimum rank r/2. Then the 
subspaces spanned by the matrices C{G\) and C(G 2 ) within 
the Schubert cell S a are at a minimum subspace distance 
dmm > r, where C(Gj), i = 1,2, are defined as above. 
Proof: The dxn matrices C{ = C(Gi), i = 1,2, have the 
matrices G\, G2, which belong to the additive group and have 
minimum rank r/2, as the d x (n — d) submatrices constituted 
of the columns indexed by b%, b 2 , • • • , b n -d- Then we have: 



d S {C u C 2 ) =2 rank 



c 2 



-2d = 2 rank 



C\ — C2 

c 2 



-2d > r. 



The last step follows from the fact that the non-zero columns 
of the difference matrix C\ — C 2 constitute an element of the 
additive group of matrices of minimum rank r/2. ■ 
Hence we can describe our two-step construction of non- 
constant dimension subspace codes with a minimum subspace 
distance d m in as follows. 

(1) Selection of Schubert cells: Select the largest subset S of 
Schubert cells subject to the following condition: For any 
S a , Sp 6 S, we have A({a}, {f3}) > d 

(2) Selection of subspaces in each Schubert cell: The sub- 
spaces within a Schubert cell are chosen such that their 
spanning matrices satisfy the following condition: The 
submatrix indexed by the characteristic columns of the 
complementary Schubert cell matrix belongs to an addi- 
tive group of minimum rank |_d m , n /2j. 

The non-constant dimension subspace code is, therefore, the 
union: [JC a , where the union is over all the Schubert cells 
S a G S. 

VI. Codes with Injection Distance as Metric 

In this section we first introduce the notion of injection dis- 
tance and discuss its applicability vis-a-vis subspace distance. 
Next we modify the construction of non-constant dimension 
subspace codes presented earlier to adapt it for the injection 
distance metric. 

A. The Adversarial Channel Model and Injection Distance 
between Subspaces: 

Recall that in the linear multicast network coding model the 
received packets at a destination are expressed as the rows of 
a matrix: Y = AX + DZ, where the matrix A is the network 
transfer matrix seen by the source packets (rows of X) and D 
is the transfer matrix seen by the error packets (rows of Z). In 
the general random network coding problem, the matrix A is 
unknown to the receiver except for a lower bound on the rank 
of the matrix. In the worst-case adversarial model investigated 
in 1 10 1, the following additional assumptions are made: 

(1) The adversary knows the matrix X and can choose the 
transfer matrix A while respecting the rank constraint; 

(2) The matrix D is arbitrarily chosen by the adversary. 



Definition 6.1: [10] The injection distance between sub- 
spaces U, V of P™ is given by: 

di(U,V) := max{Dim [/, Dim y } - Dim(E/ n V) 
= Dim(Z7 + V) - min{Dim U, Dim V}. 

The injection distance between two subspaces is thus related 
to the subspace distance between them as follows (|10|): 

di(U,V) = hi s {U,V) + hpmiU -T)imV\ (3) 

It is evident that, upto a scaling factor, the injection distance 
coincides with the subspace distance for constant-dimension 
subspace codes. It was observed in [ 1 1 that the injection 
distance gives a more precise measure of the subspace code 
performance when a single error packet can cause simultane- 
ous error and erasure. However, in a less pessimistic channel 
model than the worst-case adversarial channel, for instance, in 
the case where a single error packet can introduce either an 
error or an erasure but not both, the subspace distance is still 
appropriate. 

B. Construction of Non-constant Dimension Subspace Codes: 

First we prove a counterpart of Proposition [3] with the 
same setting and notation, for the injection distance di(U, V) 
between subspaces U, V of P™. 

Proposition 4: Let U, V be subspaces of P™ of projective 
dimensions s and t respectively; and let {fco, fci, • • • , k s }, < 
k a < n, and {Zoj^iv Jt}, < I/3 < n, be the tuples 
specifying the identity submatrices of the respective Plticker 
coordinate matrices. Then di(U, V) > [^A({k a },{lp})\, 
where A({u}, {v}) is the symmetric distance between the sets 

{«},{«}• 

Proof: From the proof of Proposition [3] we obtain: 

Dim([/fW) < (l/2)(s + t-(a + b)), where a := \{k a }\{lp}\ 
and b := \{lp} \ {k a }\. So we have: 

2d/(C7,F) := 2max{Dim [/, Dim ^ } - 2Dim([/ n V) 

> Dim U + Dim V - 2Dim(J7 n V) 

> a + b. 

As A({k a }, {lp}) = a + b, the proposition follows. 
Alternatively, from Equation [3] and Proposition [3] we have: 

2d I {U,V)>A{{k a },{lp}) + \s-t\ 

and the proposition follows. ■ 
As suggested above, to capture the effect of the injection 

distance metric in our framework, we modify the definition of 

the symmetric distance between two sets. 

Definition 6.2: For two finite sets Si , S 2 of respective 

cardinalities si, s 2 , the modified symmetric distance is defined 

as follows: 



A m (Si, S 2 ) = #{Si \ S 2 } + #{S 2 \ Si} 



si - s 2 



where #S denotes the cardinality of the set S and '| |' denotes 
the absolute value. 

Proposition 5: The modified symmetric distance is a metric. 



Proof: It is evident that A m (S , i,5 2 ) = A(5i, S 2 ) + |si - s 2 \. 
Symmetry and the fact that A m (5i,52) > with equality if 
and only if Si = S 2 follow immediately. Triangle inequality 
is satisfied because A( ) is a metric and for positive integers 
si,S2,s 3 we have: |si - s 2 | + |si - s 3 | > |s 2 - s 3 |. ■ 
Now we reformulate the two-step construction of non-constant 
dimension subspace codes with injection distance as the 
metric. From Propositions |4] and |5] the selection criterion 
of Schubert cells for the injection distance metric is read- 
ily obtained from the subspace distance case. Moreover, as 
the second step is essentially a construction of a constant- 
dimension subspace code within a Schubert cell, the transition 
from subspace distance to injection distance involves only a 
scaling factor. Hence we may state the construction of non- 
constant dimension subspace codes with a minimum injection 
distance S m i n as follows. 

(1) Selection of Schubert cells: Select the largest subset S of 
Schubert cells subject to the following condition: For any 
S a , Sp G S, A m ({a}, {/?})> 2S mm ; 

(2) Selection of subspaces in each Schubert cell: The sub- 
spaces within a Schubert cell are chosen such that their 
spanning matrices satisfy the following condition: The 
submatrix indexed by the characteristic columns of the 
complementary Schubert cell matrix is an element of an 
additive group of minimum rank 8 m i n . 

The non-constant dimension subspace code is the union: (J C a , 
over all the Schubert cells S a G S, as defined earlier. 

VII. Non-constant Dimension Ferrers-diagram 
Rank-metric Subspace Codes 

In the proposed construction of non-constant dimension 
subspace codes, if the additive group Q of d X (n — d) 
matrices is chosen to be a rank-metric code 11231 , we have a 
so-called non-constant dimension Ferrers-diagram rank-metric 
code (FDRM) subspace code. We show that the Khaleghi- 
Kschischang construction ifTTI . which has resulted in some 
of the best-known subspace codes in terms of code rate, is 
subsumed in our framework when the additive group Q is a 
rank-metric code. 

In the proposed construction, the FDRM codeword matrices 
have non-zero elements in positions which match the asterisks 
in the 'matrix' M(S a ), where S a is the Schubert cell charac- 
terized by the tuple {a} = {ai,a 2 , • • • , a<i}- First we prove a 
lower bound on the dimension of such an FDRM code, similar 
to Theorem 3 in IfTTI . in this general setting. 

Theorem 7.1: The dimension k of the largest FDRM code 
whose codewords have non-zero elements precisely in the 
positions of the asterisks in M(S a ) satisfies: k > D a — 
max{<i, n— d}(6 r — 1), where 6 r is the minimum rank distance 
of the rank-metric code and D a , the number of asterisks in 
M(S a ), is the dimension of the Schubert cell S a . 
Proof: It is known |23l that the dimension of a linear 
maximum rank-distance (MRD) code C, whose codewords are 
d x (n — d) matrices over F 9 , is given by: 

D = max{d, n — <i}(min{e?, n — d} — 5 r + 1) 



where 5 r is the minimum rank distance of C. 
Define V := w q dx{n ~ d) . The kernel of the surjective map: 
7r : V — !> V/C is a subspace of dimension d(n — d) — D. 
The largest FDRM code which matches M(S a ) is the largest 
subcode C C C subject to the constraint that all its codeword 
matrices have zeros in d(n — d) — D a fixed positions. Invoking 
the standard vector-space isomorphism between ^ q mxn and 
¥ q mn where m, n are positive integers, any codeword c G C 
may be expressed as a d(n — d) -dimensional vector over ¥ q 
with non-zero elements in D a coordinates ,«d q - Let 

(j> : V — > Yq d ^ n ^ d ^ D " be the surjective map which sets zero 
the coordinates ,«_d q . Now consider the linear map: 

$ : V — > V/C x ¥q d(n-d)-D a . v ^ (t^),^)). 

Evidently, Dim$(T/) < 2d(n — d) — D — D a . The largest 
FDRM code is the kernel of the map $, and hence, by the 
rank-nullity theorem, has dimension: 

k > d(n -d)- (2d(n — d) — D — D a ) 

= D a + max{<i, n — ei}(min{<i, n — d} — S r + 1) — d(n — d) 

= D a — max{<i, n — d}(S r — 1). 

■ 

A. Selection of Schubert Cells: 

It follows from Theorem 17.11 that the size of the FDRM 
non-constant dimension subspace codes depends directly on 
the dimension of the chosen Schubert cells. In 0111 the 
relation between the so-called profile vectors associated with 
each Schubert cell and the dimension of the FDRM codes 
constructed has been used in an algorithm for selection of the 
Schubert cells. We present a more direct approach, basing our 
choice on the dimension of each Schubert cell, as follows. 
Recall from Subsection III-CI that the dimension of a Schubert 
cell S a in the second representation characterized by the tuple 
{a} is given by: D a = nd{a) - £g a t - ( d( 2 a) ), where 
d(a) :— \{a}\ = d is the length of the tuple. Let A denote 
the set of all such distinct tuples a in lexicographic order. For 
each a G A define the choice function: 

F(S a ) =D a - max{d(a), n - d(a)}(5 r - 1). (4) 

The selection of Schubert cells proceeds via a greedy algo- 
rithm, similar to that in [11 J, which selects at each stage the 
tuple {a} with the largest F(S a ) and discards all the tuples 
j3 € A such that A({a}, {/?}) < d m ;„ for the subspace dis- 
tance metric or LiA({a},{/3})J + Li|#{a}-#{/3}|J < S rnm 
for the injection distance metric. The next stage proceeds 
by applying the same choice function to the modified set 
A := A\B, where B is the set of the discarded tuples {(3} at 
the previous stage. The selection algorithm terminates when 
A := 0. The choice function in Equation |4] which is based on 
the lower bound on the dimension of the FDRM subspace code 
given in Theorem 17.11 ensures that at every stage a Schubert 
cell which supports the largest FDRM code is chosen among 
the contending cells. The largest number of subspaces belong 
to the principal Schubert cell of the Grassmannian Q(d,n), 
with d = [n/2\ or d = fn/2]; for these values of d, the 



negative term in the choice function is also minimum for a 
given S r . Hence the algorithm proceeds by choosing either of 
these subspaces (different when n odd) in the first step. 

B. Comparison with the Khaleghi-Kschischang algorithm: 

The above algorithm and that stated in [11| differ in 
two aspects, one of which can be explained as a general 
reformulation and the other, a slight modification. 

(i) Binary Profile Vectors vs. Characteristic Tuples: In 111 II . 
the selection of Schubert cells is based on the so-called 
binary profile vectors associated with each cell. The 
Ferrers diagram representation of each profile vector then 
gives the dimension and the positions of non-zero entries 
of the appropriate rank-metric code. In our construction 
we simply use the characteristic tuple representation of 
each Schubert cell and the rank-metric code is likewise 
specified by the so-called Schubert cell "matrix" (with the 
asterisks in the "basis vectors" replacing the dots in the 
Ferrers diagrams) - a direct formulation. 

(ii) The Choice Function for the Greedy Algorithm: The 
following scoring function is used to select a new profile 
vector v in IfTTI : 

n i 

score (v, S r ) = /J/J^iVj — max{wi(u), T](v)}(d r — 1), 

»=1 3=1 

(5) 

where rj{v) — n — (wt(v) + mm ies / v \t) + 1, wt(v) the 
(Hamming) weight of the binary vector v € {0,1}™, 
v the complement of v and s(v) the support of v, 
i.e. the set whose elements denote the positions of the 
l's in v. It is evident that the sum term in the above 
scoring function computes D a , the number of asterisks 
in the 'matrix' representation of the chosen Schubert 
cell in our framework, which, in our choice function 
(Equation 0J, is the exact formula based on the char- 
acteristic tuple. Rewriting the second term in Equation 
|4] as m&x{d(a),n — (d(a) + mini {a,}) + l}(S r — 1) 
makes it equivalent to the second term of Equation [5] 
It modifies the parameter n — d(a) by subtracting the 
number of leftmost all-zero columns in M(S a ), which is 
not necessary for making a choice in our framework. 

C. Numerical Results: 

We compute lower bounds on the rates of codes constructed 
for a given minimum subspace distance (using symmetric 
distance) and for a given minimum injection distance (using 
the modified symmetric distance), from the values of the 
choice function (Equation HI for each selected Schubert 
cell. For a selected set of Schubert cells S, a lower bound 
on the code rate of a non-constant dimension FDRM 
subspace code over a finite field ¥ q is given by the formula: 

lo g 9 (Es a e5^ (5cv) ) where 

T(s a ) = f F{Sa) ' ifF ^^° 

I otherwise. 



q 


dr 


d s 


N 


LB(d s ) 


KK(d s ) 


LB(dj) 


KK(d 7 ) ) 


2 


2 


4 


9 


15.1515 


15.1518 


15.3238 


15.6245 


2 


2 


4 


10 


20.1534 


20.1534 


20.1967 


20.3294 


2 


2 


4 


12 


30.1556 


30.1557 


30.1998 


30.3346 


2 


3 


6 


13 


28.0030 


28.0030 


28.0134 


28.0263 


3 


2 


4 


7 


8.0131 


8.0145 


8.0464 


8.1331 


3 


2 


4 


8 


12.0135 


12.0135 


12.0160 


12.0311 


4 


2 


4 


7 


8.0030 


8.0032 


8.0142 


8.0522 


4 


2 


4 


8 


12.0031 


12.0031 


12.0034 


12.0068 



TABLE I 

Comparison of lower bounds on the rates of our FDRM 

SUBSPACE CODES WITH THE ACTUAL RATES OF CODES IN fTTl 



It is seen that the calculated lower bound is not simply 
a computation of the choice function in each case and then 
summing the values in the exponent. The justification for this 
modification is as follows: The choice function F(S a ) for 
a Schubert cell S a G S is zero or negative when it does 
not support a rank-metric code of the given dimensions with 
the prescribed minimum rank-distance. However, the particular 
Schubert cell still maintains the inter-cell minimum subspace 
or injection distance and so contributes at least one subspace 
codeword to the entire code. To give a simple illustration, 
consider the following construction of a non-constant dimen- 
sion FDRM subspace code in P 4 (F2) with minimum subspace 
distance 4. A set S for this code consists of Schubert cells 
specified by the tuples {1, 2}, {3, 4}, {1, 2, 3, 4, 5}, with cor- 
responding choice function values of 3, —1, —5, respectively. 
For the second tuple M(S a ) has only 2 asterisks while the last 
one has none, being the 5x5 identity matrix, obviously not 
enough to support the required rank-metric code in each case. 
However, they can still contribute a single subspace codeword 
each. 

In Table 1 we have given a few lower bounds on the code rates 
in our construction with respect to given minimum subspace 
distance (column LB(d s )) and minimum injection distance 
(column LB(d/)) for some small field sizes q and projective 
dimensions N — 1. It is seen that they are comparable with 
those of the actual codes constructed in IfTTI . which are also 
presented in the columns ' KK(d/)' (minimum injection dis- 
tance) and 'KK(d s )' (minimum subspace distance). In fact, for 
subspace distance the rates of the constructed codes actually 
coincide with the computed lower bound in half of the cases, 
while the other values show only marginal improvement. The 
difference between the rates of the constructed codes and 
the computed lower bounds is more pronounced for injection 
distance. 

VIII. Conclusion 

We have established a unified framework for the construc- 
tion of constant-dimension random network codes in terms of 
the classical description of the subspaces of a projective space 
using the Pliicker coordinates. Further we have given a general 
construction of non-constant dimension subspace codes, which 
exploits a basic representation of Schubert systems, for both 
the subspace distance and the injection distance metrics. We 
have demonstrated that our framework subsumes the construe- 



tion of ifTTl . which is corroborated by the proximity of the 
estimated lower bounds for our construction with the actual 
code rates given for small parameters in iTTTI . Some questions 
regarding the above constructions remain open: 

(i) For the non-constant dimension FDRM subspace codes, 
will an exhaustive search algorithm yield a significantly 
better code-rate than the lower bound based on a greedy 
search? Is there a systematic search procedure to achieve 
the best possible code-rate? 

(ii) Is it possible to obtain a set (or better still, a group) of 
subspaces within a selected Schubert cell larger than the 
largest Ferrers-diagram rank-metric code, which satisfies 
the minimum subspace/injection distance condition? This 
problem is a constant-dimension subspace code construc- 
tion problem with the added constraint of the Schubert 
cell structure. 

So far there have been several approaches to the decoding 
problem of constant dimension codes in general ((2]), of 
spread codes in particular, fll24l ) and employing list decoding 
methods ( 1251 . Il26ll . ll27lD . Decoding of 'punctured' non- 
constant dimension codes is discussed in 0. A possible 
decoding strategy of our non-constant dimension subspace 
codes and those in IfTTl involves two steps: 

(1) Estimating the Schubert cell to which the transmitted 
subspace belongs. In ifTTll . so-called binary asymmetric 
codes with minimum asymmetric distance greater than or 
equal to the minimum injection distance and binary codes 
with minimum Hamming distance greater than or equal 
to the minimum subspace distance have been used to esti- 
mate the profile vector corresponding to the Schubert cell 
containing the transmitted subspace. In our construction, 
a similar procedure, identifying each characteristic tuple 
with a codeword of a binary asymmetric code or that of a 
binary code with the same relations on the injection and 
subspace distances, would suffice for estimating the tuple 
corresponding to the Schubert cell. 

(2) Estimating the transmitted subset within the Schubert cell, 
which is the decoding problem for rank-metric constant- 
dimension subspace codes. 

Hence another fruitful direction of research will be to devise 
an efficient general decoding algorithm for non-constant di- 
mension projective space codes. 
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